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Abstract 

The validity of our Monte Carlo simulation procedure {the integral method) 
had been verified by the corresponding analytical procedure of which is quite 
independent of our method methodologically Also, the results obtained by 
our procedure are compared with those obtained by the different Monte Carlo 
simulation procedure {the differential method) which have been exclusively 
utilized by the different authors and the agreement between them are found 
to be well. By utilizing our Monte Carlo procedures, the validity of which 
is guaranteed in two different procedures, we investigate not only the fluc- 
tuation of high energy muons themselves but also fluctuation of the various 
quantities related to the energy losses by the muons, which are difficult to 
obtain by the differential method. Namely, we obtain fluctuation on energy 
losses of the muons, fluctuation on Cherenkov lights due to the accompanied 
cascade showers initiated by the muon and the correlations between them . 
Finally, we obtain the transition curves for Cherenkov lights in KM3 detector, 
taking into account of all possible fluctuations in the stochastic processes and 
point out the difficulty of the reliable estimation of the energy of the muons 
which are resultants of muon neutrino events in the KM3 detectors. 
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1. Introduction 



The fluctuation in high energy muon's behavior may play an important 
role in the analysis of muon neutrino events for KM3 detector deployed in 
the Antarctic, the ocean and the lake [l[ Q Q. As far as the treatment 
of the range fluctuation of high energy muons by the Monte-Carlo method is 
concerned, there exist two independent approaches. The one is the differential 
method in which the muons concerned are pursued in step by step way jij] j^] 
|!1 10 1 11]. In this method, the quantity of v cut is introduced so as 



to separate the continuous parts from the radiative parts in the stochastic 
processes in order to save the time for computation. The other is the integral 
method in which the interacti on p oints of the muons and their dissipated 



energy are directly determined [12| |22j and here, all the processes are treated 



in the stochastic manner without the introduction of v cu t- These two methods 
are independent form each other, but are logically equivalent, giving the same 
results as for the muons' behaviors (see Figures 4 to 6). However, it should be 
noticed that the energy determination of the high energy muons made by the 
measurement of the Cherenkov lights which are produced by the accompanied 
cascade showers. These cascade showers are generated from the stochastic 
processes, such as bremsstrahlung, direct electron pair production and photo 
nuclear nteraction which are initiated by muons concerned. 

2. Range fluctuation of the (ultra-) high energy muons and indi- 
vidual behavior of the muons 

2.1. The physical meaning of "no fluctuation" 

The average energy loss by high energy muon is usually described as, 

*^ = a{E ll ) + b (E„) ■ E„ (1) 

where a is the term due to ionization which is free from fluctuation and b 
is the term due to stochastic processes which may be origins of fluctuations. 
The latter is divided into three parts. Namely, 

b(E lt ) = b brems + b d . p + b nucl (Ep) , (2) 

where 6fe rems , bd. p and b nuc \ are the corresponding terms due to bremsstrahlung, 
direct electron pair production and photo nuclear interaction, respectively. 
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In the treatment of the average energy loss, each b term is defined as, 

J\J PVmax 

h (Efj) = - v 

where v max and v min are the maximum and the minimum fractional energies 
due to their kinematical limits. The physical meaning of Eq.flTJ is that the 
muons concerned disipate energy uniquely, being defined by Eq.(j2]), namely, 
fluctuations are not included in Eq.(j2]). 

In Figure [lj we give the b terms due to different processes in water. 



doj (v, E p 
dv 



■ dv, 



(3) 




In Figure [21 we give the relation between a(-E^) and b(Efj). As the b-terms 
essentially are of stochastic character, it is seen from the figure that the 
stochastic processes become effective above ~ 1 TeV. Therefore, we must 
treat the muon's behavior in stochastic manner above ~ 1 TeV. Below ~ 1 
TeV we may treat muon's behavior in the non-stochastic manner. 
Then, the range of the muon is uniquely determined by Eq.(j3J). 

f E ° d x 

R= I 3 , (4) 

where E min is the minimum energy for observation and E is the primary 
energy of the muon. Through the present paper E min denotes the minimum 
energy E min among the energies for observation (E obs ) and it is taken as 1 
GeV. Thus, R defined by Eq. (jJJ) gives the effective range of the muon without 
fluctuation. Exactly speaking, R is the muon range where the fluctuation 
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effects in the stochastic processes are neglected. The physical meaning of "no 
fluctuation" is that the muons in the stochastic processes lose their energies 
in the form of the effective energy loss defined by Eq.Q. 

In Figure [31 we give the effective range defined by Eq. (J4]) together with 
the average ranges of the muons in which the fluctuation effects are exactly 
taken into account. (see, discussion in the later sections). 

It should be noticed from the figure that the effective ranges without 
fluctuation are different from the average range of the muons in which the 
fluctuation is considered, which Lipari and Stanev[5j already pointed out. 
Really, the real average ranges are smaller than those of effective range be- 
yond one standard deviation above 10 13 eV as shown Figure [3j 

2.2. Physical quantities with fluctuation 

In the differential method, many authors divide 



all stochastic processes into two part, namely, the continuous part and radia- 
tive part in their Monte Carlo simulation in order to consider the fluctuation 
in both the range and energies of the muons and introduce v cut to save time 
for computation, while we treat all stochastic processes as exactly as possible 
without the introduction of v cut . The validity of our Monte Carlo method had 
been checked by the corresponding analytical method which is methodologi- 
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cally independent of the Monte Carlo procedure and the essential structure 
of our Monte Carlo method is described in the Appendix. Furthermore, we 
check the validity of our method, comparing our results with the correspond- 
ing result by the differential methods, which are shown in (12.2.11) . 

2.2.1. The comparison of our results with others 

Our survival probability for high energy muon is defined as, 

ry I J-, V T? \ ^through 

F (E obs ,X,E p ) = — — — — , (5) 

-t* sample y^p) 

where E p , X, and E Q b s , denote the energy of primary muon, the point for 
observation and the minimum energy of the muon at the point for observa- 
tion, respectively. N samp i e (E p ) denotes the total sampling number of muons 
and N through (E obs , X) denotes the number of muons concerned with energies 
above E obs which pass through the observation point X. 

We compare our results by the integral method with the different authors' 
results by the differential method in the following. Lipari and Stanevj^] give 
the survival probabilities as the functions of the depths for 10 11 eV to 10 18 eV 
incident muons in water and partly standard rock, the minimum energy of 
which is taken as 1 GeV. We obtain the corresponding results by the integral 
method and compare our results with Lipari and Stanev's in Figure HI Also, 
Klimushin et al[8| give the survival probabilities for primary energy of 10 13 eV 
to 3 x 10 16 eV. We obtain the corresponding results to them and compare our 
corresponding results with the results by Klimushin et aljsj in Figure |5j 



Furthermore, Kudryavtsev[lO[j gives the energy spectrum of the muon due 



to primary energy of 2 TeV at 3 km in water. We obtain the corresponding 




1000 



Figure 4: The comparison of Figure 5: The comparison Figure 6: The comparison 
our result with that of Lipari of our results with that of of our results with that of 
and Stanev Klimushin et al Kudryavtsev 
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results to him and compare our corresponding results with his results in 
Figure O The agreements between the different authors' result obtained by 
the differential method and our results obtained by the integral method are 
well as shown in Figures 3 to 5, taking into account of the slight differences 
in the cross sections utilized between the different authors' and ours. It is 
seen from these figures that the validity of our integral method is guaranteed 
by the differential method due to the different authors. 

As already explained, all the authors due to the differential method whose 
results are compared with ours divide the stochastic processes into two parts, 
namely, the radiative part and the continuous part to perform the Monte 
Carlo calculation for the study on the fluctuation of the muon behaviors. 
For the purpose, they introduce v cut by which they separate the radiation 
processes from the continuous part and they study the fluctuation effect of 
the muon in the radiative part only. 



~dE~ 




dx 


rad 



dE 

dx 

A E 



+ 



soft 

Vcut 



dE 



dx 



dv ■ v 



rad 

dv 



N F 



dv.v da{v : E) (6) 



I'cut 



dv 



Such treatment is logically correct only as far as we are interested in the 
muon behaviors, because the energy loss by the muon with single primary 
energy is exactly taken into account in their treatment irrespective of any 
v cut . However, if we are interested in Cherenkov radiation responsible for 
all stochastic processes through the accompanied cascade showers, then, the 
methods adopted by these authors are not adequate for the study on such 
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the purpose, on which we discuss in later section, (see, section H]) 

2.2.2. Survival probabilities and their differential energy spectra at different 
observable depth 

In Figures [7] to |9j we give the survival probabilities for different observable 
energies with primary energies of 10 12 , 10 15 and 10 18 eV, respectively. In Fig- 
ured we give the minimum observation energies 10 9 eV, 10 10 eV and 10 n eV, 
respectively. In Figure^ we give them 10 9 eV, 10 10 eV, 10 n eV, 10 12 eV, 10 13 eV 
and 10 14 eV, respectively. In Figure [HI we give them, 10 9 eV to 10 17 eV, re- 
spectively. Each sampling number in Figure 7 to 9 is 100,000. It is seen from 
the figures that the survival probabilities become remarkably large as their 
primary energies increase. 

In Figures [TD] to |T21 we give the differential energy spectra for primary 
energies, 10 12 , 10 15 and 10 18 eV at different depths, respectively. It is seen 
from the figures that the energy spectrum at the initial stage are of delta- 
function and they shift as the mountain-like deforming their shape in the 
intermediate stage and finally, they disappear as the results of the delta- 
function type again. Each sampling number in Figures 10 to 12 is 100,000. 




Figure 10: Energy spectrum Figure 11: Energy spectrum Figure 12: Energy spectrum 
at the different depths, initi- at the different depths, initi- at the different depths, initi- 
ated by 10 12 eV muons atcd by 10 15 eV muons ated by 10 18 eV muons 



2.2.3. Range Distribution of Muon 

All processes, such as bremsstrahlung, direct electron pair production 
and photo nuclear interaction are of stochastic ones and, therefore, one can- 
not neglect their fluctuation essentially. The muons propagate through the 
matter as the results of the competition effect among bremsstrahlung, direct 
electron pair production and photo nuclear interaction. 
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Figure 13: Range distributions for 10 12 eV to 10 18 eV muons. The minimum observation 
energies are taken as 10 9 eV. Each sampling number is 100,000. 



In Figure [T3| we give P (R; E p ), the probabilities for the range distribution 
with primary energies, 10 13 , 10 15 and 10 18 eV, respectively. It is clear from 
the figures that the width of the range distribution increases rapidly, as their 
primary energy increases. Also, as the primary energy decreases, the width 
of range distribution becomes narrower and approaches to the delta function- 
type, the limit of which denotes no fluctuation. Their average ranges and 
their standard deviations are given Table [TJ Then, the range distribution 
can be well approximated as the normal distribution in the following. 

where E p , < R > and a are primary energy, the average value of ranges and 
the standard deviations, respectively. 

In order to examine the nature of the stochastic processes in each pro- 
cess further, we compare the real range distribution in which each stochastic 
process (bremsstrahlung, direct electron pair production and photo nuclear 
interaction) are taken into account as the competition effect with the hypo- 
thetical range of the muon due to each stochastic process. 

In addition to the real range distributions, we give the hypothetical range 
distributions for 10 13 , 10 15 and 10 18 eV, in which only one cause among three 




S 



Tabic 1: The average values and values of their standard deviations of the range distri- 
butions of muons which are expressed as the normal distribution. These numerical values 
are given in both the km and g/cm 2 for primary energies from 10 eV to 10 18 eV. 



E p 
[eV] 


< R > 
[g/cm 2 ] 


< R > 
[km] 


a 

[g/cm 2 ] 


a 
[km] 


10 11 


4.75 xlO 4 


4.75X10" 1 


2.39xl0 4 


2.39x10^ 


10 12 


2.51xl0 5 


2.51x10° 


5.32xl0 4 


5.32X10" 1 


10 13 


7.18x10 s 


7.18x10° 


1.97xl0 5 


1.97x10° 


10 14 


1.23xl0 6 


1.23X10 1 


3.35xl0 5 


3.35x10° 


10 15 


1.73xl0 6 


1.73X10 1 


4.35xl0 5 


4.35x10° 


10 16 


2.22xl0 6 


2.22X10 1 


5.16xl0 5 


5.16x10° 


10 17 


2.70xl0 6 


2.70X10 1 


5.86xl0 5 


5.86x10° 
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stochastic processes is taken into account as shown in Figures [14] to [T6j Here, 
the symbol of [d. p. only] in these figures means the hypothetical range dis- 
tribution in which only direct electron pair production is taken into account 
and the bremsstrahlung and nuclear interaction are neglected. The symbols 
of [brems only] and [N.I. only] have similar meaning to that of [d. p. only]. 
From the shapes of the distribution and their maximum frequencies for dif- 
ferent stochastic processes in Figures [14] to [16], it is clear that energy loss 
in the direct electron pair production is of small fluctuation, while both the 
bremsstrahlung and photo-nuclear interaction are of bigger fluctuation and 
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the fluctuation in photo nuclear interaction becomes remarkable when com- 
pared with bremsstrahlung as primary energy increases. 

2.2.4- The diversity of individual muon behavior 

In Figures 17 to 25, we show the diversities of the muons' behaviors for the 
same primary energy of muons with regard to their ranges (or their energy 
losses). Also, in Table 2 and Table 3, we summarize the characteristics 
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Figure 17: The energy losses Figure 18: The energy losses Figure 19: The energy losses 
with the shortest range for with the average-like range with the longest range for 
10 12 eV muons. for 10 12 eV muons. 10 12 eV muons 




Figure 20: The energy losses Figure 21: The energy losses Figure 22: The energy losses 
with the shortest range for with the average-like range with the longest range for 
10 15 eV muons. for 10 15 eV muons. 10 15 eV muons. 



of the events in Figure 17 to 25. Let us examine the characteristics of the 
events in Figures 17 to 25, combined with Tables 2 and 3. The energy losses 
for given primary energy of 10 12 eV due to different stochastic processes are 
given as the function of the depths traversed. 

In Figure 17, we give the case with the shortest range. The ordinate 
denotes the ratio of the energy loss to the primary energy and D,B,N de- 
note the causes of energy losses, namely, direct electron pair production(D), 
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Figure 23: The energy losses Figure 24: The energy losses Figure 25: The energy losses 
with the shortest range for with the average-like range with the longest range for 
10 18 eV muons. for 10 18 eV muons. 10 18 eV muons 



bremsstrahlung(B) and photon nuclear interaction(N), respectively. The ab- 
scissa denotes the depths where the interactions concerned occur. The con- 
tinuous lines (red lines) represent the muon energies at the corresponding 
depths. 

It is seen from Figure 17 and Tables 2 and 3 that this muon lose ~10% 
of the primary energy at ~2 meters from the starting point and the almost 
of the primary energy by the bremsstrahlung at ~62 meters after several 
negligible experiences of energy losses due to both the direct electron pair 
production and photo nuclear interaction. 99.9% of the total energy loss is 
lost by only two bremsstrahlung. This is a good example of showing the 
catastrophic energy loss due to the bremsstrahlung. 

On the other hand, in Figure 19, we show the case with the longest range. 
The range given by Figure 19 , ~2700 meters, is far longer compared with 
~62 meters with the shortest range. It is clear from the figure and tables 
that 69.7% of the total energy is lost by 288 direct electron pair production. 
Only 29.5% of the total energy is lost 8 bremsstrahlung and the contribution 
from photo nuclear interaction is negligible. 

In Figure 18, we show the case with the average-like range. The definition 
of the average-like range denote the case whose range is the nearest to the 
average range which is obtained from the total number of the events (100,000 
sampled events). This case shows that 36.2% of the total energy is lost by 232 
direct electron pair production, 34.4% is lost by 4 photo nuclear interaction 
and 29.4% is lost by 3 bremsstrahlung, while in the real average case, 52.6% 
of the total energy is lost by 239 direct electron pair production, 33.7% by 
bremsstrahlung and 13.7% by photo nuclear interaction. 

In Figures 20 to 22, we show the similar relations for 10 15 eV muons 
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Table 2: The details of the characteristics on the muons with the shortest range, the 
aver a ge-like range ,thc longest range and the average range. 





Range 


Energy loss 


Number of 


Energy loss 


Number of 


Energy loss 


Number of 


E p = W v2 eV 


[km] 


by brems 


interaction 


by direct pair 


interaction 


by nuclear 


interaction 


A vpracp*> 


2.39x10° 


1.09xlO u 


4.67x10° 


1.70x10" 


2.39xl0 2 


4.43X10 10 


3.37x10° 


Shortest 


6.23x10-- 


9.81xlO u 


2 


4.73x10 s 


4 


3.99xl0 8 


1 


Average- like 


2.51x10° 


1.17x10" 


3 


1.44x10" 


232 


1.37x10" 


4 


Longest 


9.29x10° 


6.09xl0 10 


8 


1.44x10" 


288 


1.66x10° 


2 


E t = 10 15 eV 




<Average> 


1.72X10 1 


3.38x10" 


4.64X10 1 


4.96x10" 


6.56xl0 3 


1.61x10" 


5.31X10 1 


Shortest 


1.68x10° 


8.59x10" 


5 


8.87xl0 13 


628 


5.20x10" 


4 


Average-like 


1.73X10 1 


4.23x10" 


35 


4.75x10" 


6363 


9.67xlO Ki 


56 


Longest 


3.16X10 1 


1.82x10" 


70 


7.50x10" 


10981 


5.84xl0 13 


91 


E p = 10 18 eF 




<Average> 


3.17X10 1 


3.24 xlO 17 


1.04xl0 2 


4.59x10" 


2.48xl0 4 


2.17x10" 


1.66xl0 2 


Shortest 


8.27x10° 


7.04x10" 


33 


2.91x10" 


7285 


4.77xl0 15 


40 


Average-like 


3.19X10 1 


3.76x10" 


95 


3.51x10" 


22395 


1.64x10" 


173 


Longest 


5.48X10 1 


7.52xl0 16 


190 


7.61x10" 


45367 


2.73x10" 


299 



as shown in Figures 17 to 19 and Tables 2 to 3. The shortest range, ~1.6 
kilometer (Figure 20), is far shorter compared with the longest one, ~22 kilo- 
meters(Figure 22). It is seen from Figure 20 and Tables that bremsstrahlung 
plays a decisive role as the cause of catastrophic energy loss, too(~80% 
of total energy at ~840 meters). 85.9% of the total energy is lost by 5 
bremsstrahlung, 8.87% by 628 direct electron pair productions and 5.20% by 
4 photo nuclear interactions. In Figure 22, we give the case with the longest 
range. Here, large number of the direct electron pair production with rather 
small energy play an important role, similarly as shown in Figure 19. Here, 
75.7.8% of the total energy is lost by 10981 direct electron pair production, 
18.4% by 70 bremsstrahlung and 5.9% by 91 photo nuclear interaction. In 
Figure 21, we give the case with the average like range. Here, 47.8% of the 
total energy is lost by 6363 direct electron pair productions, 42.5% by 35 
bremsstrahlung and 9.72% by 56 photo nuclear interactions, while in the 
real averages, 49.8% of the total energy loss is lost by 6560 direct electron 
pair productions, 34.0% by 46.4 bremsstrahlung and 16.2% by 53.1 photo 
nuclear interactions. 

In Figures 23 to 25, we show the similar relations for 10 18 eV muons as 
shown in Figures 17 to 19. The case with shortest range in Figure 23 has 
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Table 3: The Ratios of energies transferred from brcmsstrahlung, direct electron pair 
production and photo nuclear interaction to the total energy loss. 
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4.22X10" 1 


3.94X10- 1 
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6.78xlQ- 2 


6.86X10- 1 
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strong contrast to that with the longest range. The manner of the energy 
loss in Figure 23 is drastic with two big catastrophic energy loss due to 
bremsstrahlung, while that in the Figure 25 is moderate with no catastrophic 
energy loss. The shortest range , ~8 kilometers, is far shorter compared with 
the longest range, ~54 kilometers. It is seen from Figure 23 and Tables that 
bremsstrahlung is a decisive role as the cause of catastrophic energy loss. 
70.4% of the total energy is lost by 33 bremsstrahlung, 29.1% by direct 
electron pair productions and 0.477% by 40 photo nuclear interactions. In 
Figure 25, we give the case with the longest range. Here, 68.8% of the total 
energy is lost by 45367 direct electron pair productions, 24.6% by 299 photo 
nuclear interactions and only 6.78% by 190 bremsstrahlung in the complete 
absence of catastrophic energy losses. In Figure 24, we give the case average- 
like range. Here, 39.88% of the total energy is lost by direct electron pair 
productions, 42.2 % by 95 bremsstrahlung and 18.4% by 173 photo nuclear 
interactions, while, in the real average's, 45.9% of the total energy is lost 
by 24800 direct electron pair productions, 32.4% by 104 bremsstrahlung and 
21.2% by 166 photo- nuclear interactions. Thus, it is concluded that the 
diversity among muon propagation with the same primary energy should be 
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noticed. 



3. Cherenkov lights production due to both the energy losses by 
the muon(naked muons) and the accompanied cascade showers 
initiated by the muons concerned 

It should be noticed that the energy losses by high energy muons are 
never measured directly. Usually, in high energy neutrino astrophysics ex- 
periments in water (ice), they are measured via Cherenkov lights which are 
produced not only by the muon itself, but also accompanied by the cascade 
showers due to bremsstrahlung, direct electron pair production and photo 
nuclear interaction, all of which are generated by the parent muons. When 
the muons traverse the matter, they lose their energies by bremsstrahlung, 
direct electron pair production and photo nuclear interaction in addition to 
ionization. These stochastic processes produce cascade showers whose pri- 
mary particle is a photon in bremsstrahlung and is an electron(a positron) in 
the direct electron pair production and photons decayed from ttq and others 
in photo nuclear interactions. These accompanied cascade showers produced 
by these stochastic processes are twisted around the traversing muons and 
these showers produce Cherenkov light. In this section, we discuss various 
quantities obtained from high energy muons, imaging the one-cubic kilome- 
ter detector for high energy neutrino astrophysics, something like Ice Cube 
in the Antarctic. 

We simulate exactly not only the both the interaction points and dissi- 
pated energies due to all stochastic processes, but also simulate exactly the 
accompanied cascade showers themselves due to these stochastic processes 
for the calculation of Cherenkov lights. Namely, the total Cherenkov lights 
due to both the muon itself and the accompanied cascade showers are ex- 
actly simulated. Here, we adopt the one- dimensional cascade showers under 



Approximation B [13[ as cascade showers. Concretely speaking, we simulate 
cascade showers as exactly as possible so that the segments of the simulated 
electrons in the cascade showers are decided in both their locations and ener- 
gies, which produce finally Cherenkov lights with the attenuation coefficient. 

3.1. The ratio of the Cherenkov lights production due to the accompanied 
cascade showers to the total Cherenkov lights 

In Figure |26l we give the ratios of the Cherenkov lights due to the ac- 
companied cascade showers to the Cherenkov lights due to (the accompa- 
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- 100 meter 

- 200 meter 

- 300 meter 



Figure 26: Ratio of Cherenkov lights due Figure 27: Observation points 
to the accompanied cascade showers to for the correlation between energy 
the total Cherenkov light losses due to muon and the pro- 

duced Cherenkov lights 



nied cascade showers plus muon itself as the functions of the depth). Near 
~ 10 11 eV, most of Cherenkov light comes from muon itself(naked muon). 
Near ~ 10 12 eV, about half of total Cherenkov light comes from muon itself 
muon. Near ~ 10 13 eV, 90% of the total Cherenkov light comes from the 
accompanied cascade showers. Above 10 14 eV, most of the total Cherenkov 
light is due to the accompanied cascade showers. 

3.2. The correlation between the total Cherenkov lights and the corresponding 
energy losses 

We examine the following correlations at observation points as shown in 
Figure [271 In Figures [28] to [35], we give the correlation diagram at obser- 
vation points, such as, 200, 500 and 1000 meters from the incident points 
in the case of 10 11 to 10 16 eV muons, respectively, between the energy losses 
due to the stochastic processes in addition to ionization and Cherenkov lights 
which are produced by both the muon themselves and their accompanied cas- 
cade showers. The attenuation coefficient is considered in the propagation 
Cherenkov lights. Here, [ energy loss ] denotes that the energy dissipated by 
the muon while traversing through some distance (for example, 200 meters) 
due to bremsstrahlung, direct electron pair production and photo nuclear 
interaction in addition to ionization loss. [Cherenkov lights] denote the mea- 
sured Cherenkov lights at some depth (for example, 200 meters) which are 
produced by not only the muon itself, but also, the accompanied cascade 
showers due to the possible stochastic processes taking into account of the 
attennation effect of Cherenkov lights. 
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Figure 30: Correlation dia- Figure 31: Correlation dia- Figure 32: Correlation dia- 
gram for 10 13 eV muons gram for 10 eV muons gram for 10 15 eV muons 



In Figure [28] we give the correlation between the energy loss and Cherenkov 
lights up to 200 meters from the starting point initiated by 10 n eV muon(naked 
muon). As seen from Figure [26] the most Cherenkov lights are produced by 
muon itself and the smaller part is produced by the accompanied cascade 
shower. Red points denote the correlation at 200 meters It is understood 
from the figure that more dense region may come from the muon itself and 
more weaker dense region may come from accompanied cascade showers. We 
cannot observe Cherenkov lights at both 500 meters and 1000 meters, be- 
cause the energy of 10 11 eV is too small to detect at 500 meters and 1000 
meters. 

In Figures [29] to [35] green points and blue pints stand for the correlation 
at 500 meters and 1000 meters, respectively. In Figure [29], we give the similar 
diagram as shown in Figure [28] and ,there, we give the correlations at 200 
meters, 500 meters and 1000 meters. As seen from Figure [26] ~ half of the 
total Cherenkov lights may be due to muon itself 's origin and ~the other 
half may be accompanied cascade showers' origin. It is clear from the figure 
that the domain for the correlation at 200 meters is larger than those at 500 
meters and 1000 meters, the meaning of which shows bigger fluctuation at 
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Figure 33: Correlation dia- Figure 34: Correlation dia- Figure 35: Correlation di- 
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200 meters compared with observation points. 

The larger part of the energy losses may be initiated by the accompanied 
cascade showers. At 10 13 eV, as seen from Figure [2H1 ~ 90 % of the total 
Cherenkov lights may be produced by the accompanied cascade showers and 
the domain for the correlation at different depths begin to overlap. Above 
10 14 eV, almost of Cherenkov lights are produced by the accompanied cascade 
showers. As the developments of the cascade showers are suffered from the 
fluctuation. It is clear from Figures EH and [32] that the boundaries of the 
domains for the correlation at 10 14 eV and 10 15 eV become obscure. 

In Figures [33] to [35], we give the correlation at different depths, say 200 
meters, 500 meters and 1000 meters for 10 16 eV muon, separately and respec- 
tively. The comparison among Figures shows that fluctuation in the energy 
losses become decrease, as the depths increase. Also, it should be noticed 
from the figures that the degree of the fluctuation in the total Cherenkov 
lights is bigger than that of the energy losses. 

The fluctuation of the total Cherenkov lights is directly related to that 
of the accompanied cascade showers which start different depth having at 
different primary energies. 

3. 3. Fluctuation in the energy loss distribution for given primary energies at 
the depths, 200 meters, 500 meters and 1000 meters. 

In Figures I3B1 to [351 we give the energy loss distribution for 10 12 eV at 200 
meters, 500 meters and 1000 meters, respectively. In Figures [3H] to and 
Figures |4~2] to 144] we give the similar distributions for 10 15 eV muons and for 
10 16 eV, respectively. In these Figures, we add the normal distributions with 
the same average values and same standard deviations which are obtained 
from the real distributions in Figures [36] to HH It is clear from these figures 
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that the normal distributions for the energy losses never express the real dis- 
tributions on the contrast to the cases of range fluctuation (see, Eq.(J7])). The 
cause of the fluctuation in the energy losses comes from the compound effect 
of the both fluctuation in the interaction points due to different stochastic 
processes and the fluctuation of the energy release from different stochastic 
processes. 





■if ir.y 


oss distribution 






\ 


Av 




AVE-1.16*1t 
SD 1 03-K 


!■ ' i 










Depth 
E p - 
lODOflO 
Stof 


=200nr 
12 eV 

Sample. 
1 0K--| 































Average E= 

' SD= 


2.79x10 2 [GeV] 
1.41>:l6 J [GeVl 










DeDih=50um 
E -10 12 eV 
OOOOr fsimnlo 
Stop 1 Okm 



























Energy Loss distribution 
AVE 5:i2"iO ? [GcV] Avcrfl 9 c 
SD 1 52xl0 ? [GeV] 
Dipl i 1000m 

E-tlofeV 
loonoo Snm^ 

Stop ' 0\rr 



n energy loss[GeV] 




energy loss[GeV] 
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3.4- Fluctuation in the total Cherenkov lights quantities for given primary 
primaries energy at the depths, 200 meters, 500 meters and 1000 meters. 
The primary muons are dissipated their energy by bremsstrahlung, direct 
electron pair production and photo nuclear interaction in addition to the 
ionization loss. Thus, these stochastic processes are the origin of the accom- 
panied cascade showers which finally produce Cherenkov lights in addition 
to the Cherenkov lights due to the muon itself (naked muon). It is clear from 
the Figures H5] to [53] that the Cherenkov lights quantities thus obtained by 
stochastic processes are widely distributed due to the complicated compound 
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fluctuation effect coming from various stochastic sources. The normal dis- 
tribution whose average values and the standard deviation are taken from 
the real distributions are given with the real distributions for the comparison 
in Figures [45] to [53] It is easily understood that such normal distributions 
never reflect upon the real situations. The correlation between the energy 
losses and Cherenkov lights are obtained from the combination of the energy 
losses as shown in Figures I3"6l to BH with the corresponding Figures [4"51 to 1531 



3. 5. Transition curves of the averaged Cherenkov lights for different primary 
energies 

In figure 54, we give the average transition curve for Cherenkov lights 
for 10 12 eV muon in the KM3 detector. The sampled number of the events 
is 500. The uncertainty bars in the figures show the range within which 68 
% of the total number is included which correspond to one a in the normal 
distribution. As shown in Figures 45 to 56, Cherenkov lights distributions are 




Figure 45: Cherenkov lights Figure 46: Cherenkov lights Figure 47: Cherenkov lights 
distribution at 200 meters for distribution at 500 meters for distribution at 1000 meters 
10 12 eV nations 10 12 cV muons for 10 12 eV muons 
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far deviated from the corresponding normal distribution and ,therefore, the 
ranges for the uncertainty too large those of the normal distributions. Also, 
gradual decreases of the slopes in the average transition curves denote that 
the energy losses up to the depths cannot be neglected compared with the 
primary energies over 1 kilometers. In Figures 55 and 57, we give the similar 
transition curves to Figure 54 for 10 13 eV and 10 15 eV, respectively. Compared 
these curves with the curve in Figure 54, it is seen that the average energy 
losses are almost independent on the behaviors of the muons concerned over 
1 kilometers. However, the uncertainties around the average values become 
large compared with that in Figure 55, which denotes bigger uncertainty as 
the increase of the primary energies. It should be noticed from these figures 
that there are the uncertainties by one order of magnitude above 10 13 eV 
muons. 
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4. Discussions and Summary 

The validity of our Monte Carlo simulation procedure (the integral method) 
has been proven to be right by two different methods. One method is an- 
alytical method (see, the Appendix) and the other is different Monte Carl 
procedure (the differential method). (see, Figures 4 to 6). 
We apply our method to get the survival probabilities for different energies 
(10 12 eV, 10 15 eV and 10 18 eV). Also, we have obtained the range distributions 
for different primary energies (10 11 eV to 10 18 eV) and we have shown they 
can be expressed as the normal distribution well (see, Table 1). We have 
given the three typical muon behaviors [with the shortest range], [with the 
average-like range] and [with the longest range] for the different primary en- 
ergies (10 12 eV, 10 15 eV and 10 18 eV) to clarify the diversity among the muons' 
ranges. (Figure 17 to 25) and have summarized their characteristics in Table 
2 and Table 3. 

In addition to the fluctuation of high energy muons, we have investigated fluc- 
tuations the total Cherenkov lights produced by muons. The total Cherenkov 
lights consists of two parts, namely, The Cherenkov lights due to the naked 
muons themselves and the Cherenkov lights through the accompanied cas- 
cade showers due to bremsstrahlung, direct electron pair production and 
photo nuclear interaction due to the muons. Above 10 14 eV, almost of the 
Cherenkov lights are produced exclusively by the accompanied cascade show- 
ers(See, Figure 26). The fluctuation of the Cherenkov lights due to the 
accompanied cascade showers come from not only the fluctuations of their 
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Figure 56: Transition curve of the averaged Figure 57: Transition curve of the averaged 
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interaction points due to the different stochastic processes, but also fluctua- 
tions of the transferred energies due to these processes due to the muons and 
the fluctuation on the development of the accompanied cascade showers. We 
have simulated all possible stochastic processes related to the production of 
the Cherenkov lights as exactly possible and have obtained the dependence 
of the average Cherenkov lights on the depths traversed (the average tran- 
sition curves for the Cherenkov lights ) (Figure 54 to 57). It is clear from 
the Figures that fluctuation of the Cherenkov lights at each depth are very 
big, the ranges for uncertainties which are of one order of the magnitude. 
From these transition curves, we can estimate primary energies within some 
allowance, which is supposed to be one factor 5 order, really big uncertainty. 
Here, it is adequate for us to mention the essential difference between the 
differential me£/iod[5],[6],[7],[8],[9] and the integral method[l2],[22] as for the 
Cherenkov lights propagations. Here, we return to the fundamental equation 
for high energy muon in order to focus on the dissipated energies by the 
muons. Again, the equation is given as, 
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In the differential method, the Monte Carlo procedure is made on the hard 
part of Eq.(8) including the integral from v cut to 1. The specified value of v cut 
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is not so sensitive for the conclusions on the muon propagations, although 
usually v cut =10~ 2 to 1CT 4 . Because, the total amount of the energy loss 
by the muon concerned is the same, irrespective of v cut in the integrals of 
Eq.(7). Therefore, the average behaviors of the high energy muons obtained 
by the differential method are same as obtained by the integral method(see, 
Figures 4 to 6). However, the shape of the energy loss spectrum obtained 
by the differential method are expected to be different from those obtained 
by the integral method, because the introduction of constant v cu t for different 
primary muon energies into the differential method may make the real energy 
loss spectrum deform Q . 

The energy loss spectrum at different depths are constructed through 
the compound processes in which the fluctuations of the interaction points 
due to the different stochastic processes and fluctuations in the transferred 
energies due to the different stochastic processes are exactly considered. The 
energy loss spectrum thus constructed are the origins of Cherenkov lights 
spectrum. Detailed speaking, the Cherenkov lights measured at some depth 
are the product by the complicated aggregate of the cascade showers whose 
starting points and whose primary energies are different due to the different 
stochastic processes. The muon energy loss spectrum should be considered 
in such the framework. 

As for the muon energy spectrum, we have another matter to be con- 
sidered related to the differential method. In the differential method, one 
introduces the constant value of v cut for the Monte carlo procedure. Here, it 
should be noticed that the primary muons have energy spectrum. Suppose 
that when E is greater than E x v cut , E may lie in the radiation part of Eq.(l) 
in the differential method. However, the E with the same energy may lie the 
soft ( continuous ) part, if their primary may be higher. In other words, 
one cannot treat the energy loss spectrum as well as the muons' behaviors 
themselves in consistent manner in the differential method with the constant 
values of v cut . 

Summarized speaking, the muon loss energy spectrum obtained by the 
differential method, probably, one cannot construct the reliable energy loss 



Exactly speaking, in the integral method, v cu t is introduced in the bremsstrahlung to 
avoid the catastrophic energy loss. However, v cut is taken as v cut <C E m i n /Eo, where Eq 
and E m i n denote the primary energy and the minimum energy for observation (1 Gev). 
Usually, we adopt 10~ 4 to 10~ 10 as v cut . Consequently, we can reproduce the differential 
cross sections for bremsstrahlung exactly above E m i n . 
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spectrum by which Cherenkov lights production spectrum are obtained. 

In Figures 54 to 57, we give the average transition curves for Cherenkov 
lights as the functions of the observation points over 1 kilometers. Each 
sampling number is 500 except for 10 15 eV (100 sampling). The uncertainty 
bars denote the range within which 68% of the total events are included, 
which correspond to the normal distribution with the average values and 
standard deviations obtained in these distributions. It is easily understood 
from these figures that the uncertainties of Cherenkov lights at each depth 
distribute over one order of magnitude except the case for 10 12 eV and it 
is almost pretty difficult for us to decide the muon energies from the corre- 
sponding Cherenkov lights. Namely, this means that it is almost difficult for 
us to construct the corresponding neutrino energies directly from the muon 
energies. 

Instead, an alternative way is suggested to decide the neutrino energy 
spectrum and muon energy spectrum simultaneously, via Chrenkov lights 
Production spectrum. 

Finally, we comment on the influence of LPM effect on Cherenkov light 
production spectrum. There are two different the LPM effect. The Influence 
of the LPM effect on the muon itself can be neglected up to 10 21 eV [14], [15]. 
The LPM effect on electron and photons [16], [17], [18], [19], [20], [21] is surely 
important factor in high energy neutrino astroparticle physics, particularly 
via electromagnetic cascade shower. The LPM effect becomes effective in the 
level of the cross sections at ~ 10 15 eV in water. However, it become effective 
at ~ 10 17 eV or more in the level of the electromagnetic cascade shower in 
water. 

In the subsequent paper, we would present the alternative way for the 
construction of the neutrino spectrum via Cherenkov lights production spec- 
trum. 
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Appendix A. 



If the differential cross section for bremsstrahlung [241] . direct pair production [25 
and photo nuclear interaction 26| are denoted by a b (E , E y ) dE^, cr^ (E , E e ) dE e 
and a n (E , E h ) dE h , respectively. 

Then, the mean free paths of the muons due to bremsstrahlung, direct elec- 
tron pair production and photo nuclear interaction are given, respectively, 

as 

\ b (E Q ) = ir -^ j (A.l) 

f Ie^ a b (E , E 7 ) dE 1 

>min 

A dp (E ) = s (A.2) 

AJE ep . a dP ^0, &ep) Cl^ep 



mi n 



An (E ) = p (A.3) 

%Jl h ™* n (Eo,E h )dE h 

Also, the resultant mean free path for there radiative processes are given as 

1 111 „, 

+ x / t— i \ + V^TT ( A - 4 ) 



\otal (Eq) \ b (E ) \ dp (Eq) \ n (Eq 

The integrations for (A.l) to (A.3) are performed over kinematically allow- 
able ranges. In (A.l), E ltmin is taken to be satisfied with such a condition 
that E 1>min / Eq sufficiently smaller than E min /E , where Eq and E min denote 
the primary energy of the muon and the minimum energy of the muon for 
observation. E min is taken as 1 Gev throughout present paper. In Figure 
A.l, we give the mean free paths for bremsstrahlung, direct electron pair 
production and photonuclear interaction are given as the function of the pri- 
mary energy. 

The most important procedures in our Monte Carlo method are only two 
procedures. 

The first one is where interaction occurs. 
For the bremsstrahlung, 

The traversed distance for the interaction is determined with the use of £, 
uniform random number between (0,1), as follows. 

At = -A 6 (Eq) logZ (A.5) 

Similarly, for direct electron pair production, 

At = -X dp (E ) log£ (A.6) 
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Figure A.l: Mean free paths for bremsstrahhmg, direct electron pair production, photo 
nuclear interaction and their total processes 



Similarly, for the photo nuclear interaction, 

At = -\ h (E )log£ (A.7) 

In our Monte Carlo simulation, the energy loss due to each stochastic pro- 
cesses are sampled by the following equations with the use of £, the uniform 
random number, between (0, 1) . 
For Bremsstrahlung 

£ = L _ (A.8) 

j E ' max <Jb (Eo, Ej) dE 1 

'min 

From (jOjl . 

E J = F b (E ,0 (A.9) 
For direct electron pair production 



IeZ . a dp(E ,E ep )dE ep 

£ = Ee ( A - 10 ) 

/ e epmax a dp E ep ) dE ep 

e Pmin 

From f lATTO]) . 

E e = F d (E ,0 (A.ll) 

For photo nuclear interaction 



Ie!! a n(Eo, E h ) dE, 



h 



£ = F m '" (A.12) 

J^a n (E ,E h )dE h 



26 



From (1X12]) . 



Eh = F h (E ,^) 



(A.13) 



For the discrimination among stochastic processes in our Monte Carlo simu- 
lation let us introduce the following equations. 



L (E ) 
Zb(E ) 



1/X b (Eq 



1/X b (E ) + 1/X dp (E ) + 1/X n (E ) 

l/X b (E ) + l/X dp (E ) 



(A.14) 



(A.15) 



1/X b (E ) + l/X dp (E ) + l/X n (E ) 
In Figure |A~2| we give a flow chart for our Monte Carlo simulation. The valid- 
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ity of our Monte Carlo procedure by the integral method should be carefully 
examined in the two methods which are independent of each other. Namely, 
The first one is the comparison of our procedure with the analytical theory 
which is explained in the present Appendix and the other is the comparison 
of our procedure (the integral method) with the different kind of Monte Carlo 
procedure (the differential method) which is mentioned in the text. 

Particularly, it is the best that the results obtained by a Monte Carlo pro- 
cedure are checked by the procedures which is methodologically independent 
of the Monte Carlo procedure, reaching the same results. 

In Figure A. 3, the average energies of the muons are given as the func- 
tion of the depths under the preposition of muon energy spectrum at sea level 
with indices 2, 3 and 4, obtained by the integral method [12, 22] and they 
are compared with results obtained by the analytical theory based on the 
Nishimura-Kamata formalism in the cascade shower theory and the agree- 
ment between them [23] are very well, which surely guarantee the validity of 
our integral method. 




Figure A.3: 
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